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ABSTRACT 

A fraction of gamma-ray bursts exhibit distinct spectral features in their 
prompt emission below few 10s of keV that exceed simple extrapolations of the 
low-energy power-law portion of the Band spectral model. This is also true for 
the prompt optical emission observed in several bursts. Through Monte Carlo 
simulations, we model such low-energy spectral excess components as hadronic 
cascade emission initiated by photomeson interactions of ultra-high-energy pro- 
tons accelerated within GRB outflows. Synchrotron radiation from the cascad- 
ing, secondary electron-positron pairs can naturally reproduce the observed soft 
spectra in the X-ray band, and in some cases the optical spectra as well. These 
components can be directly related to the higher energy radiation at GeV ener- 
gies due to the hadronic cascades. Depending on the spectral shape, the total 
energy in protons is required to be comparable to or appreciably larger than the 
observed total photon energy. In particular, we apply our model to the excess 
X-ray and GeV emission of GRB 090902B, and the bright optical emission of 
the "naked-eye" GRB 0803 19B. Besides the hard GeV components detected by 
Fermi, such X-ray or optical spectral excesses are further potential signatures of 
ultra-high-energy cosmic ray production in gamma-ray bursts. 

Subject headings: cosmic rays — gamma rays: bursts — gamma rays: theory - 
radiation mechanisms: nonthermal 
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Introduction 



Gamma-ray bursts (GRBs) are potential sources of ul tra-high-energy cosmic rays (UHE- 
CRs, IWaxman E99I Ivietri Il995i iMilgrom fc Usov Il996h . While the required loc al UHECR 



emiss i yity is £ p dN v I d s p ~ 0.8 xlO 44 erg Mpc 3 yr 1 at proton energy e p ~ 10 19 eV fjWaxman Sz Bahcall 



1998: 


1 ■' — r 
Dermer 


(IDaiene et al. 



2006 



2007 ), the local GRB rate may be in the range of 0.05 — 1 Gpc 



Le fc Dermer 1 120071 ; iGuetta fc Piran 1 120071 ) . Assuming p p 



1 

2 for 

the power-law index of accelerated protons, the necessary total isotropic-equivalent energy 
in protons would be E p ~ 2 x 10 54 — 3 x 10 55 erg when integrated over e p ~ 10 9 — 10 20 eV. 
Compared to the isotropic-equivalent energy released as gamma-rays, typically E 1 ~ 10 53 
erg, this indicates that the accelerated protons must dominate the energy budget of GRBs 
in order for them to be viable sources of UHECRs. 

Although Ep/E^ > 10 — 100 in GRBs may appear physically demanding , a number 
of considerations provide some justification f lTotani I Il998t lAsano et al. 1 120091 ). First, in 
the popular internal shock model where the prompt MeV gamma-rays are attrib uted to 
synchrotron radiation from electrons accelerated in shocks within the GRB outflow (IPiran 
20051 ; iMeszaros 1 120061 ). the efficiency with which the electrons are energized may be limited 
in comparison with the total available energy, which is the dissipated fraction of the bulk 
kinetic energy conveyed predominantly by protons. Since simple Coulomb collisions are 
ineffective in transferring the energy of the protons to the electrons within a dynamical 
timescale, it has been generally assumed that this can occur through some plasma instabilities 
with efficiencies e e ~ 0.1-0.5. However, this is by no means physically guaranteed, so we 
are motivated to consider the possibility of proton-dominated GRBs with e c <C 0.1 and 
explore its consequences. Note that a high proton-to-elec tron ratio is not only ob served in 



Galactic cosmic rays and inferred in supernova remnants jAharonian et al~ 20061). but als o 



theoretically expected, at least for nonrelativistic shocks (IBlandford 



1994 



Levinson Ill996l ). 



The prompt emission spectra of GRBs are known to be generally well described by 
the so-called Band model, consisting of a hard, low-energy power-law part, a softer, high 
energ y power-law part, and a spectral peak in between in the sub-MeV range (IBand et al. 
19931 ). However, the Large Area Telescope (LAT) onboard the Fermi Gamma- Ray Space 
Telescope (Fermi) recently detected an additional, hard spectral component above ^0 .1 



GeV in the prompt phase of the short GRB 090510 flAbdo et al. 112009b 



Ackermann et al. 



20101 ) . Notwithstanding alternative models (IGranot Il2010l ; iDermer Il2010f ). if the GRB accel- 
erated ultra-high-energy protons with isotropic-equivalent luminosity L p > 10 55 erg/s, syn- 
chrotron and inverse Compton (IC) emission from an electron-posi tron pair cascade triggered 
by photopion interactions of the protons with low-energy photons ( Bottcher fc Dermer 19981; 
Gupta fc Zhang 112007 ; Asano fc Inoue 12007 ) can account for this GeV emission ( Asano et al. 
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2009^). 



Interestingly, in this same burst, the Fermi Gamma-ray Burst Monitor (GBM) observed 
a further, soft excess feature below ~ 20 keV, which appears to lie on a continuation of 
the hard, GeV power-law. A similar, even clearer X-ray excess, as well as a GeV exces s 
was also reported for the prompt emission of the long GRB 090902B (lAbdo et al. 112009a ). 
Evidenc e of extra low-energ y components was also reported in about ~ 15 % of BATSE 



bursts (IPreece et al. I Il996l ). In the optical band, several GRBs have exhibited prompt 



optical fluxes t hat are brighter than expected from simple extrapolations of their low-energy 
Band spectra f lYost et al. 1 120071 ; iPanaitescu 1120081 ). One of the most impre ssive cases was 



the e xtremely luminous optical emission of the "naked-eye" GRB 0803 19B f lRacusin et al. 
20081 ). 



Despite alternative expla nations such as an early onset of the afterglow flGhisellini et al. 



2009 



2009 



2010 



Pe'er et al. 



Kumar &: Duran 1120091) or upscattering of external/photospheric photons (IToma et al. 



20101 ). the emission from hadronic pair cascades can also potentially 



account for the excess components in the X-ray or even optical band, owing to the generally 
very wide distribution in energy of the secondary pairs and their resultant radiation. In 
addition to high-energy gamma-rays, prompt X-ray and optical emission may thus prove to 
be valuable observational signatures of UHECR acceleration in GRBs. Here we demonstrate 
this through Monte Carlo spectral modeling, focusing on the two remarkable cases of GRB 
090902B and GRB 080319B. 



Model and Methods 



Our Monte Carlo code self-consistently calculates the photon and neutrino spectra cor- 
responding to individual pulses in the GRB prompt e mission in a one- zone approximation, 
including all relevant leptonic and hadronic processes (lAsano et al. Il2009bl ). More detailed 



descriptions of the code can be found in a series of previous pu blications (jAsano 1 12005 



Asano fc Nagataki 1120061 ; lAsano Sz Inoue 1120071 : lAsano et al. 1120091 ). In this work, we choose 
not to explicitly model the MeV-range Band spectral component, whose true origin is cur- 
rently under debate and may or may not be related to electrons accelerated in internal 
shocks. Even in the case that photospheric emission or some other mechanism is more rel- 



evant , it is plausible th at coexisting internal shocks (IMeszaros fc Rees 



2010 ; Pe'er et al. 2010 ) or magnetic reconnection regions ( Giannios 201ol ) accelerate pro- 



2000: Toma et al. 



tons to ultra-high energies within the GRB outflow. Here we simply assume that a radiation 
field consistent with the observed Band spectrum is present within the emission region, which 
also contains magnetic fields with energy density Ub, and is moving out with bulk Lorentz 
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factor T at radius R from the central engine. Protons are injected in this environment with 
isotropic-equivalent luminosity L p in the form of a power-law energy distribution with a 
high-energy cutoff, oc e~ 2 exp (— e p /e Ptraax ), above e Pjmin = 10 GeV in the comoving frame. 
The maximum proton energy £ Pimax is determined by the balance between the acceleration 
timescale t acc = ^R^/c, where i?L is the proton gyroradius, and either the comoving expan- 
sion timescale t exp = R/ (cT) or the proton cooling timescale t coo i due to synchrotron and 
photomeson production losses. We choose £ = 1 as required for UHECR production b y 



GRBs, and observationally inferred for supernova remnant shocks (jUchiyama et al. 1 120071 ). 
Of the free model parameters R, T, Ub and L p , the latter two can be quantified relative to the 
photon energy density U 1 = L 1 /(4ircR 2 T 2 ) of the Band component with isotropic-equivalent 
luminosity L 7 as Ub/U^ and U p /U^ where U p = L p /(47rcR 2 T 2 ). 

The physical processes taken into account are 1) photon emission and energy losses by 
synchrotron radiation and Compton scattering including the Klein-Nishina regime for elec- 
trons/positrons, protons, pions and muons, 2) synchrotron self-absorption for electrons/positrons, 
3) 77 pair production, 4) photomeson production for protons and neutrons, 5) Bethe-Heitler 
pair production (p + 7 — >■ p + e + + e~), and 6) decays of pions and muons. 

For the results below, we do not include the effects of intergalactic 77 abso rption due to 



theextragalactic background light, since the latest observational constraints flAbdo et al 



20101 . and references therein) indicate that it should not be severe at the energies and redshifts 



of our interest here. 



3. Results 
3.1. GRB 090902B 

GRB 090902B at z — 1. 822 was one of the most energetic long GRBs with isotropic 



energy E- lso = 3.6 x 10 54 erg (lAbdo et al. Il2009d ). The detection by Fermi of an 11 GeV 



photon together with variability on timescales ~ 53 ms during the prompt phase implies a 
minimum bulk Lorentz factor r min ~ 1000. Between T + 4.6 s and T + 9.6 s where T is the 
Fermi GBM trigger time, a hard power-law component with photon index T 7 = —1.94 was 
found above 0.1 GeV, in addition to a component well fit by a Band function with low-energy 
photon index a = 0.07, high-energy photon index f3 = —3.9 and peak energy -E pea k = 908 
keV. This Band component manifests an atypically narrow energy distribution, as apparent 
in Figure [TJ Also observed was a soft excess feature below 50 keV that is consistent with an 
extrapolation of the GeV power-law down to these energies. The gray shaded area in Figure 
[1] represents the l-a confidence region of the best-fit, unfolded spectrum for GRB 090902B, 
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based on a likelihood analysis of the combined G BM and LAT data under the assumption 
of a Band function plus an additional power-law ( lAbdo et al. 112009a ). 



ef(£) [erg/cm7s] R=10 14 cm, r=1300, U n /U,=3, iyu v =l 
1 1 1 1 1 1 1 1| — 1 1 1 1 1 1 1 1| — 1 1 iiimi — 1 1 iiimi — 1 1 iiinii 
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Fig. 1. — Model spectrum for parameters listed at the top as thick red curve compared with 
observations of GRB 090902B, for which the gray shaded area represents the 1-a confidence 
region of the best-fit, unfolded spectrum for the joint Fermi GBM and LAT data. The best- 
fit Band component is shown separately as the solid black curve. Individual contributions 
of synchrotron and inverse Compton from secondary electron-positron pairs as well as muon 
synchrotron are denoted by thin blue curves as labelled, not including the effects of 77 
absorption or synchrotron self-absorption. 

Overlayed in Figure[T]is our model of pair cascade emission induced by ultra-high-energy 
protons for the parameters R = 10 14 cm, T = 1300, U p /U~, = 3, and Ub/U^ = 1. The overall 
agreement is good, except for some deviation above ~ GeV. Considering the low photon 
statistics at the highest energies for the LAT, and the fact that the reference spectrum is an 
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unfolded one wherein a simple functional form was assumed, we do not consider this to be a 
serious discrepancy. Because the GeV component is characterized by photon index ~ — 2 with 
relatively low power compared to the Band component, it can be reproduced reasonably well 
by pair cascade em ission that is dominated by synchrotron processes in a relatively strong 
magnetic field (e.g. ICoppi Ill992l ). The necessary nonthermal proton luminosity is then not 
excessive and only comparable to the Band component luminosity. These inferences are in 
notable contrast to the case of GRB 090510, whose GeV emission was harder with photon 
index ~ — 1.6 and more luminous, warranting a stronger inverse Compton contribution from 
the secondary pairs and consequently a weaker magnetic field. This entailed a lower £p. max 
and hen ce lower photopion p roduction efficiency, leading to the requirement that L p > 10 55 
erg s" 1 flAsano et al. Il2009bl ). 

The cascade emission also extends to low energies and naturally explains the X-ray 
excess. Note that a proton synchrotron model would have difficulty accounting for such 
spectra with roughly equal power across a very broad energy range. On the other hand, for 
the present parameters, we do not expect strong emission down to the optical band, as the 
spectrum turns over at < 100 eV due to synchrotron self- absorption. Moreover, the narrow 
energy distribution of the Band component and the corresponding 77 opacity precludes the 
injection of electrons/positrons with energies low enough to radiate in the optical (see thin 
blue curve in Figure [1] where self-absorption is neglected). We also mention that the high- 
energy cutoff due to 77 absorption expected at > 10 GeV is a diagnostic feature to be tested 
in future observations of similar bursts. 



3.2. GRB 080319B 



GRB 080319 at z — 0.937 was also an energetic long GRB with isotropic gamma-ray 
energy E- mo = 1.3 x 10 54 erg. Most remarkably, it was accompanied by an extremely bright, 
prompt optical burst reaching 5.3 apparent magnitude, earning the moniker "naked-eye 
GRB" . Unlike most previous prompt optical detectio ns that could be reasonably interpreted 



as emission from an ext ernal reverse or forward shock (IFox fc Meszaros 1120061 ; iRoming et al. 



20061 ; lYost et al. 1 120071 ). the unprecedentedly well-sampled light curve of this burst revealed 
sever al peaks over the first ~ 50 s that were correlated with the variabili ty in MeV gamma- 



rays ( Racusin et al. 2008 ; Stamatikos et al. 2009 ; Beskin et al. Ilioiol ). before switching 



over to a behavior more typical of afterglows. Furthermore, the early optical fluxes lie 
far above low-energy extrapolations of the concurrent, time- resolved Band spectra, clearly 



point ing to a distinct spectral component associated with the prompt phase fjRacusin et al 
20081 ). Among different possibilities, an interpretation attributing the optical and MeV 
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emission respectively to the synchrotron and s ynchrotron-self- Compton (SSC) processes have 



often been invoked (e.g. iRacusin et al. II2008I ). 



For comparison with our hadronic cascade model, we focus on the time interval between 
T + 12 s and T + 22s where T is the Swift/BAT trigger time, during which the Band 
component reaches a maximum luminosity of 1.0 x 10 53 erg s _1 . The Band parameters we 
adopt are a = —0.45, (3 = —3.5, and E pea ^ = 748 keV, noting that (3 is only constrained by 
the data to be j3 < —3.412. Figure [2] illustrates our results for the parameters R = 10 16 cm, 
r = 1000, U p /U 7 = 45, and U B /U 1 = 3. 

Synchrotron emission from the pair cascade softens the spectrum below ~ 100 keV so 
that the low-energy slope of the Band component becomes concordant with the observed 
value of a = —0.816. The cascade emission continues down into the optical band with 
some curvature and accounts well for the observed optical intensity as long as synchrotron 
self-absorption does not set in, which necessitates T > 1000 and R > 10 16 cm. Such large 
values for T and R imply a relatively low comoving photon density and thus low photopion 
production efficiency, which scales as oc R~ 1 T~ 2 . This in turn calls for a proton luminosity 
L p ~ 10 55 ergs -1 , strongly dominating that in the Band component. Note that a sub- 
stantially larger energy budget than in MeV gamma-rays alone is also unavoidable in the 
SSC model for GRB 080319B due to luminous, second-order IC emission in the GeV band 



(IRacusin et al. 1 120081 ) . 



Occurring before the launch of Fermi, the GeV properties of GRB 0803 19B remains 
largely unknown. Upper limits at th e level of ~ 10 -5 erg cm -2 s -1 above 10 GeV were 
obtained by MILAGRO (lAune 1 12010| ). This does not contradict the high-energy compo- 
nent expected in our model, consisting of synchrotron emission from secondary pairs with 
luminosity ~ 0.1L 7 and photon index ~ —2, extending up to a cutoff ~ 100 GeV due 
to internal 77-absorption (Figure [2]). The spectral shape is conspicuously different from 
the sharply peaked on e expected from second-order IC emission in the SSC interpretation 



(IRacusin et al. 1120081 ). providing an important distinguishing feature for UHECR-induced 
emission components. The overall similarity of our model spectrum for GRB 0809 13B with 
that of GRB 090902B also encourages us to search for optical signatures in future Fermi 
GRBs. 



4. Conclusions and Discussion 



Emission from pair cascades initiated by photopion production of ultra-high-energy 
protons in GRB outflows can yield spectra with roughly equal power over a broad energy 
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Fig. 2. — Model spectrum for parameters listed at the top as thick red curve compared with 
observations of GRB 0803 19B, for which the gray shaded area represents the spectrum mea- 
sured between T + 12 s and To + 22s by Swift/BAT and Konus-Wind. The contemporaneous 
optical flux observed by "Pi of the Sky" is the black diamond. The best-fit Band component 
is shown separately as the thin black curve. Individual contributions of synchrotron and 
inverse Compton from secondary electron-positron pairs, as well as muon synchrotron and 
proton synchrotron are denoted by thin blue curves as labelled, not including the effects of 
77 absorption or synchrotron self-absorption. 



range. While GeV-Te V signatures du e to hadronic processes have been discussed previously 



by many authors (e.g. iDermer II2010I ). the relevance of associated X-ray or optical features 



had not received much attention. Here we showed that synchrotron emission from hadronic 
cascades can reproduce the excess X-ray and GeV components in GRB 090902B, as well as 
the bright optical emission in GRB 080319B. Unlike the case of GRB 090510, the necessary 
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proton luminosity for GRB 090902B is only comparable to the observed photon luminos- 
ity. On the other hand, GRB 0803 19B calls for an isotropic-equivalent proton luminosity 



10 55 erg s 1 , which may appear extreme but is in fact consistent with the energetics 



requirements for GRBs to be the origin of UHECRs fjAsano et al. 1120091 ) 



The detailed time variability properties of the low-energy excess components are not yet 
known observationally, except for the well- documented optical lig ht curve of GRB 080319B 
( Racusin et al. 20081 Stamatikos et al. 2009 ; Beskin et al. 2010 ) (see also Vestrand et al. 



2004 



Blake et al. 



2004 ) . Although an in-depth discussion is beyond the scope of this pa- 



per, we may expect the variability of hadronic cascade emission to be influenced by the 
timescale of photopion production and thus qualitatively different from purely leptonic pro- 
cesses. There should also be a close connection between the excess components at low and 
high energies. The X-ray excess may be better characterized in the near future by the 
Joint Astrophysics Nascent Universe Satellite (JANUS), which will conduct prompt GRB 
observations including the 1-20 keV band. The c ontinuing development of wide-field and/or 
rapidly-slewing telescopes (IMundell et al. 1 120101 ) should bring forth a clearer picture in the 
optical band for more bursts. The broadband variability of hadronic emission components 
will be discussed at greater length in subsequent work. 

In order to explain the complex behavior of th e after glow of GRB 080319B, a two- 
component jet model was discussed in lRacusin et al. I (120081 ). with a narrower and faster jet 
dominating the early phase of the afterglow, and a wider and slower jet describing the later 
phase. Such composite models would loosen the constraints that we obtained here for our 
one-zone model. If the wider jet is the main site of proton acceleration, and if it is dense 
enough to cause p — p collisions in the prompt phase, there may be other parameter sets 
that can reproduce the optical emission. Furthermore, the jet collimation-corrected energy 
in protons can be diminished in such a composite model. 

One of the problems for the internal shock synchrotron model of the prompt emission 
is the difficulty in accounting for the observed low-energy spectra of the Band component, 
which are generally much harder than naive expectations with a = —1.5, including the 
GR Bs dealt with in this pape r. Many authors have proposed modifications to the model 



(e.g. lAsano fc Terasawa 1120091, and references there in) or alternatives such as those based on 
photospheric emission (e.g. iMeszaros fc Rees 1120001 ). Although here we have tacitly assumed 
the presence of the Band component as is observed, in the future we intend to model all spec- 
tral components self-consistently within the context of physical processes that can reconcile 
the low-energy index observations, in order to achieve a more comprehensive understanding 
of emission mechanisms together with UHECR production in GRBs. 
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